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1.0 INTRODUCTION 


This Users Guide explains how to use the new version of the 
Naval Postgraduate School Random Number Generator Package [1], called 
LLRANDOMII. It is an outgrowth of the old LLRANDOM package [1] which 
has been the subject of several modifications and enhancements during 
the period 1978 to 1981, and which was based on the work of Lewis, 
Goodman & Miller [2]. 


This document does not go into the details of how the different 
methods operate internally, but concentrates on helping the reader to 
be able to use the package. The details are discussed in a separate 
report [3]; other references are given there. However, the philosophy 
behind the package is to provide arrays of uniform and non-uniform 
random numbers which have: 

(i) known and documented statistical properties [6]: 

(ii) are efficiently computed; and 

(iii) are reliably computed. 

With point (iii) in mind, nothing in LLRANDOM was changed unless it 
significantly enhanced the capabilities of the package. 


The entire package is single precision oriented, so all the 


parameters are either full word integers or short precision floating 
point . Beside integer and floating point pseudo-random deviates the 


package can generate Normal, Exponential, Gamma, Cauchy, Poisson, and 


Geometric random deviates, or the order statistics associated with 


these random deviates. 


This LLRANDOMII package is mostly written in IBM/360 assembly 
language but uses two subroutines written in standard IBM FORTRAN. 

Since IBM/360 Assembly language is upward compatible, the program will 
run on IBM/370 machines and on the IBM 3033. Testing was done on an 
IBM/360/65 and on an IBM 3033. All the programs are independent modules 
which follow all IBM standard linkage conventions and therefore can be 
called from all IBM high level languages which use the standard linkage 
conventions (FORTRAN, PL/I, etc.). 
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2.0 DIFFERENCES FROM THE PREVIOUS PACKAGE 

The new version of LLRANOOM (LLRANDOMII) is not compatible with 
the old version. Therefore users who may want to convert to the new 
package should change their subroutine name and argument list since 
both have changed from LLRANDOM to LLRANDOMII (subroutine names In 
LLRANDOMII are all different from corresponding subroutine neunes In 
LLRANDOM). 

The new package has a new parameter (except on Integer uniform 
generators), the fifth one on the argument list, which must be Included. 
This parameter must be either Integer 0 to Indicate that the array of 
random numbers Is not to be sorted or Integer 1 to request that the 
array of random numbers be sorted In ascending order. Table 1 contains 
all the subroutine neunes for all generators available with the argument 
list Included. 

The 4th argument In all generators Indicates the multiplier to be 
used In the congruentlal generator. The table of multipliers currently 
has two choices (1 = 16807, 2 = 397204094). There are two versions of 
each generator, the second providing deviates based on a self-shuffled 
version of the basic random number generator [3]. 

Using the long multiplier results In "better" random numbers at 
the expense of an Increase In execution time (up to twice as much In 
an IBM/360. No difference In time found In the IBM 3033. The shuffling 
feature also Improves the quality of the "uniform” random numbers gen¬ 
erated, and hopefully Improves the quality of the resulting non-uniform 
numbers generated when that Is the case. Shuffling adds only a very 
small eumount of time when used and allows the generation of sequences 
of uniforms of very long cycle. Given the availability of the long 
multiplier and the shuffling feature It may have been wise to drop the 
short multiplier. However this multiplier, originally proposed by 
Lewis, Goodman & Miller [2], Is widely used and Is only marginally de¬ 
fective. Thus It can be used to rerun old simulations or to obtain In 
this package two different pseudo-random number strezuns. 

The new package does not require the Initial call to the OVFLOV' 
subroutine to establish the Interrupt environment. This Is handled 
Internally now by the LLRANDOMII programs, and represents a significant 
Improvement. 

New generators have been added for Geometric and Poisson deviates 
enhancing the Normal, Gamma, Exponential and Cauchy capabilities of 
LLRANDOM and Its extensions [5]. A new faster algorithm Is used to 
generate Gzunma deviates when a > 1. 
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3.0 REIATIONSHIP TO IMSL (EDITION 8) , CHAPTER G 

The routines In the Intemetionel Mathematics and Statistics 
Libraries (IMSL) Edition 8 (designed by P.A.W. Lewis & J. Gentle) 
Implement substantially the same algorithms as are Implemented in 
LLRANDOMII. However the implementation in the IMSL Library is meant 
first of all to be portable across the computing systems which IMSL 
supports, and secondly is oriented toward generation of random numbers 
one at a time. Thus for efficiency the LLRANDOMII pac)cage should be 
used. 


Timings for the LLRANDOMII routines are given in Table 2 below 
for the IBM 3033. Some timings are given for IMSL subroutines. In 
all cases timings per deviate are based on generation of arrays of 
10,000 deviates. Timings for LLRANDOM were given in [1]. LLRANDOMII 
on an IBM 3033 is roughly eight times as fast as LLRANDOMII on the 
IBM 360/67. 

TABLE 2 

TIMING COMPARISONS AND SUBROUTINE EQUIVALENCES 



LRANDOMII 



IMSL 


VARIATE 

IBM 3033 



IBM 3033 


INTEGER UNIFORM 

LINT 

(2) 


GGUD (16) 


REAL UNIFORM 

LRND 

(3) 


GGUBS (9) 


NORMAL 

LNORM 

(7) 


GGNML (44) 


EXPONENTIAL 

LEXPN 

(7) 


GGEXN (23) 




.8 

(54) 

ra* .8 

(91) 

GAMMA(CHI-SQUARE) 

LGAMA ( 



GGAMR ( 



^0-2 

.5 

(18) 

^0-2.5 

(105) 


fP“ 

.8 

(5) 

fP-.8 

(21) 

GEOMETRIC 

LGEOM i 



GGEOT { 



'P-. 

98 

(16) 

'‘P-.98 

(21) 



8 

(5) 

fP-.8 

(29) 

POISSON 

LPOIS ( 



GGPON ( 



^P=20. 

(7) 

^p-20. 

(323) 

CAUCHY 

LCCHY 


(7) 

GGCAY 

(34) 


• Time in microseconds (10 secs) per deviate based on generation 
of arrays of 10,000 deviates at a time is aiven in parenthesis. 
The Assembly Language LLRANDOMII routines are usually 3 to eight 
times as fast as the FORTRAN IMSL subroutines. 

• For times to generate deviates one at a time add approximately 

7 microseconds to the values for LRANDOMII subroutines, 12 to 20 
microseconds to IMSL, subroutines. 

• Use of the shuffled generators adds about 5% to the times for 

LRANDOMII subroutines. , . 
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4.0 USING THE PACKAGE 

For most^of the uses Initialization of the "seed" IX is the 
only thing required in addition to the call to the appropriate sub> 
routine. There are some cases though when the user may want to 
generate several streams of numbers and then he must be aware of 
some rules to follow if he is concerned about being able to repeat 
the exact sequence of numbers later on. The "seed" is the starting 
value for the (integer) uniform random number generator and it must 
be initialized to any value between 1 and 2’^-l. It is updated 
automatically on every call to the generators so that next time the 
subroutine is called the sequence of numbers produced will continue 
at the point it ended last time . Normally, the user should never 
change the seed in the course of the program. If the user wants dif¬ 
ferent streeuns of numbers at different points in his program, then he 
should use different seed naunes and initialize them with different values 
for the different calling points. 


In the new version of LLRANDOM all the generators are independent 
of each other since each generator is now an independent module with 
its own separate storage. However, using the same generator at different 
points in the same program will still cause some interaction in the 
case of "shuffled" generators. This is because only one copy of the 
generator is loaded by the link-editor and the shuffling table there¬ 
fore becomes the same for the different calling points. Then, in this 
case, if repeatability is important the user must make sure that the 
"seeds" and the number of deviates requested at every calling point 
remain the same. Section 4.1 shows an example of a case where the 
sequences do not match in different runs of the program due to a "small" 
modification made to it. As Indicated before this happens only when 
using shuffled generators. 

4.1 EXAMPLES 

a) To generate 10,000 geometric deviates with parameter p>0.80, 
using the short multiplier, no shuffling, no sorting. The 
Geometric generator produces deviates X according to the 
following form of the Geometric distribution: 













v«r(x)-p/(l-p)^. 


f (k)-P{X-k}-p’'(l-p), 
which hai moamntm 
B(X)-p/(l-p), 

The FORTRAN program may look Ilka thia 
REAL *4 A(IOOOO), P 
Intager IX,N 
IX-15987 
N-10000 
P-.80 


Call LGEOM (IX,A,N,1,0,P) 

End. 

Thus there is no sorting and the Multiplier 1 is used. 

b) Generate two batches of ten samples each of Poisson deviates. 
Batch No. 1 with samples of sise Nl«5000 each and Batch No. 2 
with samples of size M2-3500 numbers. Both batches are 
from the Poisson distribution with X > 3.7, using multiplier 
2, shuffled and no sorting. 

The Poisson distribution has the form 

f(k)-P{x-k}-xW(Jc:) 
with mcxnents 

E(x)*X, var(x)«X. 
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The FORTRAN program may be eet up like thlei 

Real *4 A1 (5000) ,A2 (4000), AL 

Integer IXl, 1X2, Nl, N2 

IXl-14872 

1X2-932789 

Nl-5000 

N2-3500 

AL-3.7 

• 

e 

e 

Do 100 ISAMPL-1,10 

Calling sequence for siunples of Batch No. 1 
CALL SPOIS (IX1,A1,N1,2,0,AL) 

Calling sequence for samples of Batch No. 2 
CALL SPOIS (IX2,A2,N2,2,0,AL) 

CONTINUE 

END 

This example shows one case where the same shuffled 
generator Is used in different points In the same program, 
and Interaction between the two calls Is created. If the 
same sequences of numbers have to be repeated In a future 
run of the program, the user must make sure that 1X1,1X2, 

Nl and N2 remain the same. For example. If the program is 
run again with N2-4000 the shuffling table (which is shared 
by the two calls) will be updated more often now in the 
second call causing the numbers from the first call to be 
different. 

As indicated before this happens only when the same sub¬ 
routine name is used in both calls and they are of the 
shuffled type. 
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Ai far as the no-ahuffle type generator! being able to 
regenerate the tame aequencea of number a, all that is 
required ia to uae the aame value for the "aeed" in all 
different runa of the program. Different neunea for the 
aeed variable are atill neceaaary to avoid interaction 
between the different calla, in the caae of multiple 
calla in the aame uaer program. 

4.2 NOTE FOR PL/1 USERS 

When uaing the LLRANDOMll package from PL/1 programa the ap¬ 
propriate aubroutinea being uaed muat be declared aa external entries 
with data typea Binary Fixed (31) and Binary Float (21) for the 
parametera to be paaaed that correapond, respectively, to INTEGER 
and REAL *4 in Fortran. For example, if the LGAMA generator is used 
in a PL/1 program, the following declaration must be included in the 
program: 

DCL LGAMA EXTERNAL ENTRY (BIN FIXED(31) , (*) BIN FLOAT(21), 

BIN FIXED(31), BIN FIXED(31), BIN FIXED(31), BIN FLOAT(21)) 

OPTIONS(ASSEMBLER INTER); 


5.0 AVAILABILITY 

The LLRANDOMII package is distributed in 9 track 1600BPI magnetic 
tape containing two partitioned data sets. The first contains the 
source IBM Assembly and Fortran code and the second contains object 
code that can be used to construct a TXTLIB under VM/CMS or as a 
library in the link-edit step when used in batch mode. 

To obtain a copy nend a scratch tape to 


Naval Postgraduate School 
Monterey, CA 93940 
Operations Research Dept. 
Attn: Professor P.A.H. Lewis 
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